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Bayesian approach to the first excited nucleon state in lattice QCD * 
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We present preliminary results from the first attempt to reconstruct the spectral function in the nucleon and A 
channels from lattice QCD data using the maximum entropy method (MEM). An advantage of the MEM analysis 
is to enable us to access information of the excited state spectrum. Performing simulations on two lattice volumes, 
we confirm the large finite size effect on the first excited nucleon state in the lighter quark mass region. 



In quantum field theory, the spectral func- 
tions (SPFs) of the two-point correlation func- 
tion is expected to expose rich physical informa- 
tion which is not only for the ground state, but 
also for excited states. Nevertheless, the recon- 
struction of SPF, A(uj) from given Monte Carlo 
data of the Euclidian time correlator: C?(t) = 
/ dijjA(ijS) exp(— cji) is a typical ill-posed problem. 
The maximum entropy method (MEM) is a useful 
method to circumvent such ill-posed problem by 
making a statistical inference of the most proba- 
ble image of A{lij) based on Bayesian statistics. 

Recently the MEM analysis is widely employed 
on various problems in lattice simulations after 
the first success in our research area [|l[D. As 
for the light hadron spectroscopy, the CP-PACS 
collaboration analyzed their own high-precision 
quenched lattice QCD data using the MEM. They 
show the reliability of the MEM through checking 
consistency between the standard analysis and 
the MEM analysis after the continuum extrapola- 
tion . However above applications were carried 
out only for mesonic hadrons. In this article, we 
apply the MEM analysis to lattice QCD data for 
both spin-1/2 and spin-3/2 baryons in order to 
study the excited state spectrum of baryons. 

We are interested in a long standing puzzle re- 
garding the excited state spectrum of the nucleon, 
namely, the level order of the positive-parity ex- 
cited nucleon 7V'(1440) and the negative-parity 
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nucleon A^*(1535). The pattern of the level or- 
der between positive and negative-parity excited 
states can be found universally in the A, S and 
flavor-octet A channels. Recent quenched lattice 
QCD studies show that the wrong ordering be- 
tween N' and N* actually happens in the rela- 
tively heavy-quark mass region Thus, we ad- 
dress a serious question whether or not the level 
switching between N' and N* would occur in the 
light- quark mass region. Finding of this possi- 
bility might be directly connected to the under- 
standing of the mysterious Roper resonance. 

We remark that the simulation for the light- 
quark mass requires large lattice size since the 
"wave function" of the bound state enlarges as 
the quark mass decreases. Once the "wave func- 
tion" is squeezed due to the small volume, the ki- 
netic energy of internal quarks raises and thus the 
total energy of the bound state should be pushed 
up. This is an intuitive picture for the finite size 
effect on the mass spectrum. Such effect is ex- 
pected to become serious for the (radial) excited 
state rather than the ground state. Indeed, ref.[|| 
reported that the N' mass in the light-quark re- 
gion is significantly heavier than the mass extrap- 
olated from the heavy quark region. Their lattice 
simulation was performed on relatively small vol- 
ume {L « 1.5fm). 

To study the finite size effect, numerical simu- 
lations are performed on two lattice sizes 16^ x 32 
and 24^ x 32. We generate quenched QCD con- 
figurations with the standard single-plaquette ac- 
tion at /3 = 6.0 (fl-i « 1.9GeV). The quark prop- 
agators are computed using the Wilson fermion 
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Figure 1. Spectral functions in the A^(l/2+) 
channel and the A (3/2+) channel. 

action at four values of the hopping parameter k, 
which cover the range M-^jMp « 0.69 — 0.92. Our 
preliminary results are analyzed on 352 configura- 
tions for the smaller lattice (L w 1.5fm) and 300 
configurations for the larger lattice (L « 2.2fm). 

We use the conventional interpolating op- 
erators, eabc{uJCj5db)uc for the nucleon and 
£abc{'u-JCjf^Ub)uc for the A respectively. Corre- 
lators constructed from those operators are sup- 
posed to receive contributions from both positive 
and negative-parity states. More details of the 
parity projection are described in ref.0|. 

For the A correlator G'^^, we need the spin pro- 
jection to extract desired spin state since there are 
contributions from both J = 3/2 and 1/2 states. 
At zero spatial momentum, the A correlator with 
spatial Lorentz indices is expressed in the form 
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where contributions from negative-parity states 
are omitted for simplicity. Defining = 
jTT{'jijjGjl), one can obtain each amplitude; 
Gy^ and G^^^ from appropriate combinations of 
Gij . In this article, we only analyze the = 
1/2+ part of the N correlator and the = 3/2+ 
part of the A correlator. 

Wc define the dimensionless SPF for baryons 
through A{uj) = p{ijj)uj^ In the contin- 

uum perturbative QCD, the asymptotic (renor- 
malized) value of for the large renormal- 

ization scale IGeV) can be evaluated within 
the one-loop level 
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Here the value of ci = 5/128 (1/10) and C2 = 
71/12 (52/45) for our used N (A) operator in 
the MS-scheme can be retrieved from the number 
quoted in ref. In our MEM analysis, the prior 
knowledge (default model; mo) is deduced from 
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Zg^" denotes 
the renormalization factor of the specified baryon 
interpolating operator on the lattice. For the N 
operator, the tadpole improved perturbative re- 
sult is give by (1 — 0.73Q!y)/8 in the chiral limit 
We roughly choose mg = 4 x 10~^ for N and 
8 X 10~^ for A according to the above procedure. 
The results are insensitive to the variation of mp 
around this choice. 

The MEM analysis is performed by using data 
up to the half of the temporal extent except for 
the source location. Most of the results presented 
here use N^) ~ 0(600) and WmaxO ~ 27r. The 
detail procedure of the MEM can be found in ||^ . 

In Fig.l, we show SPFs in the nucleon and A 
channels on the larger volume (L w 2.2fm). The 
solid, dashed and dotted lines are for k=0.155, 
0.153 and 0.1515 respectively. In both channels, 
there are one sharp peak, one broad peak and 
two bumps at each k. Two bumps at each n have 
large overlap with the ones at different k, while 
the peak positions are relatively shifted to the 
left as K, increases toward the chiral limit. Those 
states might be the unphysical bound states of a 
physical quark and two doublers, which have been 
found in the mesonic case The additive sim- 
ulations at different lattice spacing are necessary 
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Figure 2. Study of the finite size effect in the 
nucleon channel at k = 0.155 and 0.1515. 

to confirm this speculation. 

In order to examine the finite size effect, we 
compare SPFs of the nucleon on the larger volume 
with the ones on the smaller volume at k = 0.155 
and 0.1515 in Fig. 2. The crosses on each peak 
or bump represent the statistical significance of 
SPF obtained by the MEM. We see a large finite 
volume effect on the second peak for the lighter 
quark (k — 0.155) as compared with the first peak 
and two bumps. It indicates that the (physical) 
excited state is significantly affected by the finite 
size effect in comparison to the ground state. 

We plot the N and N' masses, which corre- 
spond to the peak positions of first two peaks, as 
a function of the pion mass squared in Fig. 3. The 
errors are estimated by the jackknife method. We 
mention that the N masses are quite consistent 
with the ones determined from the single expo- 
nential fits. In addition, for the N', the heaviest 
two points are also good agreement with the cor- 
responding results of ref . Q| . Taking a simple lin- 
ear extrapolation in Fig. 3, we find AfAr=:0.51(l) 
and Mjv'=0. 90(10) in lattice units. Our results 
may be compared with the previously published 
results for the N* at the same lattice spacing; 
Mn- = 0.85(5) D for lattice size L « 1.5fm with 
domain wall fcrmions and 0.89(2) § for lattice 
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Figure 3. The ground state (x) and the first 
excited state (o) masses versus the pion mass 
squared in the nucleon channel. 

size L w 2.2fm with clover fermions. We find that 
the level spacing between N* and N' reduces sig- 
nificantly in the chiral limit. However the level 
switching between them might not happen in lat- 
tice simulations with L ^ 2.2fm. 

We have applied the maximum entropy method 
to lattice QCD data for both spin-1/2 and spin- 
3/2 baryons to study the positive-parity excited 
state spectrum. We succeeded in extracting SPFs 
for baryons as well as mesons. Based on the sys- 
tematic analysis utilizing two lattice sizes, we con- 
firmed the large finite size effect on the first ex- 
cited nucleon state in the light quark mass region 
originally pointed out in ref.Q. 
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